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Abstract 

Recently, Smajic et al. [1] published an article on numerical struc- 
tural optimizations of two-dimensional photonic crystals using two differ- 
ent classes of optimization algorithms i.e, deterministic for local searches 
and stochastic for global. In this comment we reexamine some of their 
conclusions regarding the global stochastic search strategies. It is con- 
cluded that the steps taken to increase the efficiency of the optimization 
of the test problem chosen by Smajic et al, which was selected to keep 
the search space tractable, can be misleading when applied to intractable 
problems. 

1 Introduction 

First of all, when evaluating the performance of an optimization strategy the 
class of optimization problems at hand play a very crucial part. In 1997 Wolpert 
and Macready [2] published the 'No free lunch theorems for optimization'. The 
theorems proof that for any algorithm, any elevated performance over one class 
of problem is exactly paid for in bad performance over another class. In other 
words, it does not exist no better nor worse optimization algorithm with respect 
to its average performance on all possible classes of problems. In Ref. [1], Smajic 
et al. analyze the performance of four different global search strategies by 
optimizing a power divider coded by 12 binary parameters. The 12 parameters 
give rise to a total of 4096 possible configuration. The relatively small size of the 
problem was chosen to keep the search space tractable, and by doing so, they 
limit themself to a specific class of problems. If this is not done carefully this 
class might be inadequately represented in a global perspective and the method 
will lack generalization. 
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2 An example 



As a counter example the problem of optimizing a doped PC lens has been 
chosen. The number of lattice sites has been increased to 182, which results 
in a considerable increase in the size of the search space. Now a problem with 
2 182 = 6.1 x 10 54 possible configurations has to be processed. Because of the 
immense number of solutions it is impossible to do a brute-force evaluation of 
all these configurations making the optimization process indispensable. Please 
refer to Ref. [3] for any further details about this specific problem. To design 
the lens structure a design tool based on a multiple scattering theory [4], as a 
direct solver, and a genetic algorithm (GA) [5] for global optimization, was im- 
plemented. This approach is similar to one of the four global search algorithms 
analyzed in section 4 of Ref. [1]. The authors of Ref. [1] was "disillusioned" 
by the poor performance of the standard GA for solving a 12 bit-parameters 
problem of designing a power divider implementing the GA with a population 
of 45 individuals. As a countermeasure they used a micro-GA that works with 
multiple runs of a standard GA implemented with a much smaller population. 
In this way they improved the performance of the optimization process consid- 
erably. Here we have repeated the some of the steps taken by Smajic et al. and 
shown that there conclusions lack the generality of the concept " Optimization 
of photonic crystal structures" . 

Figure n picture the performance of different GAs implemented to optimize 
the 182 bit-parameters lens problem. The GAs are implemented with different 
population sizes. Figure^, and b correspond to the so called micro-GA with 
population sizes of 5 and 20, respectively, while Fig. ^ an( A d correspond to 
normal GA optimizations with populations of 100 and 200 respectively. Each 
optimization was carried out five times to get an idea of their stochastic behavior 
and each run was stopped after 50 '000 evaluations. It is clearly demonstrated 
that for this particular problem, a so called micro-GA (figCJi and b) performs 
bad and a standard GA with a large population (fig^ 3 an d c) is considerable 
more robust and effective. These results contradict the conclusions drawn in 
Ref. [1] from the 12-bit power divider design problem. 

3 Conclusions 

In summary, the optimization or inverse-design of photonic crystal structures 
is an challenging problem. To find a global search algorithm for solving the 
wide class of problems consisting in effective design of PC devices is a difficult 
task that has to be tackled with rigor. To analyze the performance in a generic 
way of any proposed scheme of optimization it is important to choose a general 
problem or even to solve a set of problems. In our opinion, this has not been 
done in Ref. [1] and, therfore, their conclusions can be misleading. 
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List of Figure Captions 

Fig. 1. GA popultion dcpcndecy. The best fitness in the population at every 
generation is printed out for four different GAs. The four GAs differ in popu- 
lation sizes and are applied to a problem coded with 182 bit-parameters. The 
following population sized were used 5, 20, 100 and 200 in a), b), c) and d) re- 
spectively. Each GA is executed 5 times corresponding to the five lines in each 
graph. The average fitness value after 50'000 evaluations and for 5 independent 
runs are for a) 3.34, b) 3.69, c) 4.06 and for d) 4.25 
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Figure 1: GA popultion dependecy. The best fitness in the population at ev- 
ery generation is printed out for four different GAs. The four GAs differ in 
population sizes and are applied to a problem coded with 182 bit-parameters. 
The following population sized were used 5, 20, 100 and 200 in a), b), c) and d) 
respectively. Each GA is executed 5 times corresponding to the five lines in each 
graph. The average fitness value after 50'000 evaluations and for 5 independent 
runs are for a) 3.34, b) 3.69, c) 4.06 and for d) 4.25 
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